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Abstract 

The strong coupling constants among members of the heavy spin-3/2 baryons 
containing single heavy quark with light pseudoscalar mesons are calculated in the 
framework of the light cone QCD sum rules. Using symmetry arguments, some struc- 
ture independent relations among different correlation functions are obtained. It is 
shown that all possible transitions can be described in terms of one universal invariant 
function whose explicit expression is Lorenz structure dependent. 
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1 Introduction 



The heavy baryons have been at the focus of much attention both theoretically and exper- 
imentally during the last decade. The heavy quarks inside these baryons provide windows 
helping us see somewhat further under the skin of the nonperturbative QCD as compared to 
the light baryons. Since a part of polarization of heavy quark is transferred to the baryon, so 
investigation of polarization effects of these baryons can give information about the heavy 
quark spin. Moreover, these baryons provide a possibility to study the predictions of heavy 
quark effective theory (HQET). Besides the spectroscopy of heavy baryons, which have 
been discussed widely in the literature, their electromagnetic, weak and strong decays are 
very promising tools to get knowledge on their internal structure. In this decade, essential 
experimental results have been obtained in the spectroscopy of heavy baryons. The 
[ \~] antitriplet states A+, 5+ 5° [ A+(2593), S+(2790), S0(2790)] as well as the |+ [ 
1^] sextet states fig, Sg, [fi*, S*, S*] have been observed [1]. Among the S-wave bottom 
baryons, only the A^, E^, E^, and 0,^ have been discovered. It is expected that the 
LHC will open new horizons in the discovery of the excited bottom baryon sates [2] and 
provide possibility to study the electromagnetic properties of heavy baryons as well as their 
weak and strong transitions. The experimental progress in this area stimulates intensive 
theoretical studies (for a review see for instance [3,4] and references therein). Theoretical 
calculations of parameters characterizing the decay of the heavy baryons will help us better 
understand the experimental results. 

The strong coupling constants are the main ingredients for strong decays of heavy 
baryons. These couplings occur in a low energy scale far from the asymptotic freedom 
region, where the strong coupling constant between quarks and gluons is large and pertur- 
bation theory is invalid. Therefore, to calculate such coupling constants a nonperturbative 
approach is needed. One of the most reliable and attractive nonperturbative methods is 
QCD sum rules [5]. This method is based on QCD Lagrangian and does not contain any 
model dependent parameter. The light cone QCD sum rules (LCSR) method [6] is an ex- 
tended version of the traditional QCD sum rules in which the operator product expansion 
(OPE) is carried over twists rather than the dimensions of the operators as in the case 
of traditional sum rules. In the present work, we calculate the strong coupling constants 
among sextet of the heavy spin-3/2 baryons containing single heavy quark with the light 
pseudoscalar mesons in the framework of the LCSR. Using symmetry arguments, we show 
that all allowed strong transitions among members of these baryons in the presence of light 
pseudoscalar mesons can be expressed in terms of only one universal invariant function. 
Note that the coupling constants of heavy spin-1/2 baryons with pseudoscalar and vec- 
tor mesons have been recently calculated in [7,8]. The heavy spin 3/2-heavy spin 1/2 
baryon-pseudoscalar meson and heavy spin 3/2-heavy spin 1/2 baryon-vector meson cou- 
pling constants have also been calculated in the same framework in [9,10]. It should be 
mentioned here that the couplings of the heavy baryons with mesons is first calculated in 
[11] within the framework of HQET. 

Rest of the article is organized as follows. In section 2, we derive some structure in- 
dependent relations among the corresponding correlation functions, and demonstrate how 
the considered coupling constants can be calculated in terms of only one universal function. 
In this section, we also derive the LCSR for the heavy spin-3/2 baryon-light pseudoscalar 
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meson coupling constants using the distribution amplitudes (DA's) of the pseudoscalar 
mesons. Section 3 is devoted to the numerical analysis of the related coupling constants 
and discussion. 

2 Light cone QCD sum rules for the coupling con- 
stants of pseudoscalar mesons with heavy spin— 3/2 
baryons 

In this section, the LCSR for the coupling constants among heavy spin-3/2 baryons and 
light pseudoscalar mesons are derived. We start our discussion by considering the following 
correlation function: 

Il^. = tj d'xe^P' {V{q) \T {r],{x)nM}\ 0) , (1) 

where V^q) is the pseudoscalar-meson with momentum g, 77^ is the interpolating current 
for the heavy spin-3/2 baryons and T denotes the time ordering operator. The correlation 
function in Eq. (1) is calculated in two different ways: 

• in terms of hadronic parameters called the physical or phenomenological representa- 
tion, 

• in terms of QCD degrees of freedom by the help of OPE called theoretical or QCD 
representation. 

Matching then these two representations of the same correlation function, we obtain 
the QCD sum rules for strong coupling constants. To suppress contributions of the higher 
states and continuum, we apply the Borel transformation with respect to the momentum 
squared of the initial and final states to both sides of the sum rules. 

We start our calculations by considering the physical side. Inserting complete sets of 
hadrons with the same quantum numbers as the interpolating currents and isolating the 
ground states, we obtain 

^ ^_ (0M0)|i?2(p))(i?2(p)P(g)|i?i(p + g))(i?i(p + g)|r;.(0)|0) ^ 

{p^ — ml) [{p + g)2 — ml] 

where \Bi{p + q)) and \B2{p)) are the initial and final spin-3/2 states, and mi and m2 are 
their masses, respectively. The dots in Eq. (2) represent contributions of the higher states 
and continuum. It follows from Eq. (2) that, in order to calculate the phenomenological 
part of the correlation function, the following matrix elements are needed: 

{Bi{p + q)\r],{0)\0) = Xb.Mp + Q) , 
{B2{p)V{q) \Bi{p + q)) = QB^B^vU^iph^u^ip + q) , (3) 

where A^^ and are the residues of the initial and final spin-3/2 heavy baryons, gBiB2V 
is the strong coupling constant of pseudoscalar mesons with heavy spin-3/2 baryons and 
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Uf^ is the the Rarita-Schwinger spinor. Using the above matrix elements and performing 
summation over the spins of the Rarita-Schwinger spinors defined as 

in principle, one can find the final expression of the correlation function in phenomenolog- 
ical side. However, the following two principal problems are unavoidable: 1) all Lorentz 
structures are not independent; 2) not only the heavy spin-3/2, but also the heavy spin-1/2 
states contribute to the physical side, i.e., the matrix element of the current r]^, sandwiched 
between the vacuum and the heavy spin-1/2 states, is nonzero and determined in the 
following way: 

(0 \ri^\ B{p, s = 1/2)) = AiAp^ - mj^p, s = 1/2) , (5) 

where the condition 7/^?]^ = is imposed. To remove the contribution of the unwanted heavy 
spin-1/2 baryons and obtain only independent structures, we order the Dirac matrices in a 
specific way and eliminate the ones that receive contributions from spin-1/2 states. Here, 
we choose the ^fii)i^ulb ordering of the Dirac matrices and obtain the final expression 

Tlau — 1—^ ^\r/^^^^^fo^ QvT {qi^vMi^ + othcr structurcs with 7„ at the beginning and 

\p — ni\)\yp ^ qY — m\)\\ 

7i,75 at the end, or terms that are proportional to (p + q)y or p^ , (6) 

and to calculate the strong coupling constant QbiBqV, we choose the structure Qf^uMl^^ 
which is free of the unwanted heavy spin-1/2 states. 

Now, we proceed to calculate the correlation function from QCD side. For this aim, 
we need to know the explicit expression for the interpolating current of the heavy spin- 
3/2 baryons. In constructing the interpolating current for these baryons, we use the fact 
that the interpolating current for this case should be symmetric with respect to the light 
quarks. Using this condition, the interpolating current for the heavy baryons with J — 3/2 
containing single heavy quark is written as 

ri, = Ae^''^[{qtC^,q',)Q^ + {q',Cj,Q')q^, + (Q''C^,ql)q',} , (7) 

where A is the normalization factor, and a, b and c are the color indices. In Table 1, we 
present the values of A and light quark content of heavy spin-3/2 baryons. 

Before obtaining the explicit form of the correlation functions on the QCD side, we first 
try to obtain relations among the correlation functions of the different transitions using 
some symmetry arguments. Next, we show that all possible transitions can be described 
in terms of only one universal invariant function. We follow the approach given in [7-10], 
where the coupling constants of heavy spin-1/2 baryons with light mesons as well as spin 
3/2 baryons-spin 1/2 baryons- light mesons vertices are calculated (see also [12-16] for the 
couplings among light baryons and light mesons). Here, we should stress that the relations 
which are presented below are independent of the choice of Lorentz structures. We start our 
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Table 1: The light quark content qi and q2 for the heavy baryons with spin-3/2 

discussion by considering the — )> S^^'tt^ transition. The invariant function describing 
this strong transition can be written in the following general form: 

n^6°^s*07rO ^ ^^o_^ni(it, d, h) + g^0M^^{u, d, h) + g^^uTl^iu, d, b) , (8) 

where Qj^Ouu, 9n0dd ^'^'^ 9w%b show coupling of the meson to the uu, dd and bb states, 
respectively. The interpolating current of tt" meson is written as 



(9) 



u.d 



where g'TrOuu = —QnOdd = Q-w^b — 0- The invariant functions Hi, II'^ and 112 describe the 
radiation of 7r° meson from d and b quarks of heavy baryon, respectively, and they 
can formally be defined as: 



ni(M,d,6) = {uu\Lft.f\Q) , 

n'i(M,rf,6) = (dd|sfs^°|o> , 

Mu^d^b) = (66|S*°S*°|0> . 



(10) 



Prom the interpolating current of baryon, we see that it is symmetric under the exchange 
•<->■ d, so n^(ii, 0?, b) — ni(ci, M, b) and we immediately obtain 
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7i 



ni(M, d, b) — ni((i, M, 6) 



where under the SU2{2)f limit, ]J^b°^'^b°'^° = 0. Now, we proceed to obtain the invariant 
function responsible for other transitions containing the 7r° meson. The invariant function 
for E^"*" — >■ Hl'^n^ transition can be obtained from the — >■ E^°7r° case by making the 

replacement d ^ u, and using the fact that rj/j,'' — v2r]n'' , from which we get. 



4ni(M,'u,6) = 2('um|S*+E*+|0> . 



(12) 



The factor 4 on the left hand side appears due to the fact that each S^"*" contains two u 
quark, hence there are 4 possible ways for radiating 7r° from the u quark. From Eq. (11), 
we get 



V2Ui{u,u,b). 



(13) 
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Similar arguments lead to the following result for the — > 7r° transition: 

n^r^s*-7r° _ -^Ui{d, d, b). (14) 

Consider the strong transition El~^^^ — >■ S^~^°V°. The invariant function for this decay 
can be obtained from the E^° — >■ E^°7r° case using the fact that rj]^'' — 7]^^ {d ^ s) and 
rfj? — ViJ.'' (u ^ s). As a result, we get 

nH*-->S*-7rO ^ Lni((i,s,6) . (15) 

v2 

Now, we proceed to find relations among the invariant functions involving charged tt^ 
mesons. We start by considering the matrix element (^dd [E^'^E^"! O), where d quarks from 
the and E^*^ form the final dd state, and u and b quarks are the spectators. The matrix 
clement (^ud |E^^E^'^| O) explains the case where d quark from E^*^ and u quark from E^"*" 
form the ud state and the remaining u and b are being again the spectators. From this 
observations, one expects that these matrix elements be proportional to each other. Our 
calculations support this expectation. Hence, 

^nf-^nl+n- ^ (^rf|s*+S*°|0> = \/2((irf|EfEf |0> = y2ni((i,M,6) . (16) 

Making the exchange u -ir^ d in Eq. (16), we obtain 

^ |E*-S*°| 0) = y/2{uu |E*°S*°| 0> = V2Ui{u,d,b) . (17) 

Similarly, one can easily show that 

nE*+->E-.+ ^^ni(ci,M,6), 
nsr->E*o.- ^ y^u,iu,d,b) , 

= Ui{d,S,b) , 

jjEi-^Ef.- ^ ni(M,s,6) . (18) 



Calculation of the coupling constants of the members of heavy spin-3/2 baryons to other 
pseudoscalar mesons can be done in a similar way as we did for the meson. Here, we 
shall say that in our calculations, we ignore the mixing between t] and rj mesons and only 
consider rjg instead of physical rj meson. The interpolating current for rjs meson has the 
following form 

Jm = [ui5U + d-i^d - 2s75s] , (19) 



we see that 
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For instance, consider the — >■ S^°?78 transition. Following the same lines of calculations 
as in the 7r° meson case, we immediately obtain, 

ns^o->sf^8^ 1 [ni(M,rf,6) + ni(rf,M,6)] . (21) 

The invariant function responsible for the — > S°*r78 transition can be written as: 

U-b ^=6 ^8 = g^^^^Ui{u, s, h) + g.r,^ss^i{u, s, h) + g^^if)l2{u, s, h) 

= -^[niKs,6)-2n;Ks,6)] 

= -^[niKs,6)-2ni(s,M,6)] . (22) 
For the remaining transitions containing the r^g meson we obtain 

V6 
v6 

^ l=[u,(d,s,b)-2u,(s,d,b)] , 
n^r^^^''^ = -^ni(s,s,6) . (23) 

v6 

We also find the following relations for transitions involving K mesons: 

n=?°-^^»^° = n^r'-H-^» = ni(d,t.,6), 
n^r^E-K- ^ n^-^Hri^+ = ni(i.,d,6) , 

^U-b--^l-''°^V2U,{s,s,b), 



(24) 



The expressions for the charmed baryons can easily be obtained by making the replacement 
6 — > c and adding to charge of each baryon a positive unit charge. 

So far we have obtained all possible strong transitions among the heavy spin-3/2 baryons 
with pscudoscalar mesons that arc described in terms of only one invariant function IIi. This 
function can be calculated in deep Euclidean region, where —p'^ — )> +00 and — {p+qY — > +c>o 
using the OPE in terms of DA's of the pscudoscalar mesons as well as light and heavy quark 
propagators. In obtaining the expression of Hi in QCD side, the nonlocal matrix elements of 
types {P{q) \q{x)rq{0) \ 0) and {P{q) \q{x)Gn^q{0) \ 0) appear, where F is any arbitrary Dirac 
matrix. Up to twist-4 accuracy, these matrix elements are determined in terms of the DA's 
of the pscudoscalar mesons. These matrix elements as well as the explicit expressions of 
DA's are given in [17-19]. 
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In calculation of the invariant function Hi, we also need to know the expressions of 
the light and heavy quark propagators. The light quark propagator, in the presence of an 
external gluon field, is calculated in [20]: 



Sg{x) 



du 



2 12 

16'k'^x'^ 
I In 



G^^{ux)a^ 



-x^A^ 



ux^Gf,^{ux)Y 



(25) 



where 7^ ~ 0.577 is the Euler constant, and A is the scale parameter. In further numerical 
calculations, we choose it as A = (0.5 -j- 1) GeV (see [21, 22]). The heavy quark propagator 
in an external gluon field is given as: 



du 



[rrir 



;G''''{ux)a^, + 



u 



,(26) 



where Sq^^{x) is the free heavy quark operator in coordinate space and it is given by 



Q 



[X) 



■x^ 



47r2 



— I 



x^ 



K2{mQV-x'^) , 



(27) 



where Ki and K2 are the modified Bessel function of the second kind. 

Using the explicit expressions of the heavy and light quark propagators and definitions of 
DA's for the pseudoscalar mesons, we calculate the correlation function from the QCD side. 
Equating the coefficients of the structure gij,uMl5 i^om both sides of the correlation function 
and applying the Borel transformation with respect to the variables and {p + qY in order 
to suppress the contributions of the higher states and continuum, we get the following sum 
rules for the strong coupling constants of the pseudoscalar mesons with heavy spin-3/2 
baryons: 



9B1B2V — 



^Bi 



Hi 



(28) 



where and M| are the Borel mass parameters correspond to the initial and final heavy 
baryons, respectively. The contributions of the higher states and continuum are obtained 
by invoking the duality condition, which means that above the thresholds Si and S2 the 
double spectral density p'*(si, S2) coincides with the spectral density derived from QCD side 
of the correlation function. The procedure for obtaining double spectral density from QCD 
side and subtraction of higher states and continuum contributions are explained in detail 
in [23], which we use in the present work. 

The masses of the initial and final baryons are equal to each other, so we take Mf = 
M| = 2M2 and the residues A^i and are calculated in [24]. The explicit expression for 
Hi function is quite lengthy and we do not present its explicit form here. 
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3 Numerical results 



This section is devoted to the numerical analysis of the sum rules for strong coupling 
constants of pseudoscalar mesons with heavy spin-3/2 baryons. The main input parameters 
of these LCSRs are DA's of the pseudoscalar mesons which are given in [17-19]. Some 
other input parameters entering to the sum rules are {qq (2 GeV)) = —(274^^7 MeVY and 
H^{2GeV) = Uml/imu + m^) = (2.43 ± 0.42) GeV [25], {ss) = 0.8{uu), (0 | ^a.G^ | g) = 
(0.012 ± 0.004) GeV\ ml = (0.8 ± 0.2) GeV^ [26], U = 0.131 GeV, fx = O.^IG GeV and 
/,3=0.13 GeV [17]. 
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Figure 1: The dependence of the strong coupling constant for the — >■ S^°7r° transition 
at several fixed values of sq. 

The sum rules for the strong coupling constants include also two auxiliary parameters: 
Borel mass parameter and continuum threshold sq. These are not physical quanti- 
ties, hence the result for coupling constants should be independent of them. Therefore, we 
should look for working regions of these parameters, where coupling constants remain ap- 
proximately unchanged. The upper limit of is obtained requiring that the contribution 
of the higher states and continuum is small and constitutes only few percent of the total 
dispersion integral. The lower bound of is obtained demanding that the series of the 
light cone expansion with increasing twist should be convergent. These conditions lead to 
the working region 15 GeV^ < < 30 GeV"^ for the bottom heavy spin-3/2 baryons and 
4 GeV^ < < 10 GeV^ for the charmed cases. In these intervals, the twist-4 contribu- 
tions does not exceed (4-6)% of the total result. Our analysis also shows that contribution 
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Bottom Baryons 




Charmed Baryons 




41±9 (4.0 ± 0.5) 




22±6 (5.3 ± 0.6) 


g^b ^^b ^ 


75±15 (8.0 ±0.7) 




42±10 (11.0 ±0.7) 


g-b ^^b ^ 


65±14 (7.8 ± 0.8) 




37±8 (10.2 ± 0.6) 




73±16 (8.0 ±0.8) 




37±9 (9.8 ± 0.8) 




45±10 (5.8 ± 0.6) 




25±6 (7.8 ± 0.6) 


g^t°^^fvs 


19±4 (2.7 ±0.3) 


g^^^^^c^ris 


11.5±1.8 (3.2 ±0.4) 




95±22 (11.5 ± 1.8) 




o2±10 (13. G ± 0.7) 



Table 2: The values of the strong coupling constants for the transitions among the heavy 
spin-3/2 baryons with pseudoscalar mesons. 

of the higher states and continuum is less than 25%. The continuum threshold Sq is not 
totally arbitrary but it is correlated to the energy of the first excited state with the same 
quantum numbers as the interpolating current. Our calculations show that in the interval 
(mfi + 0.4)^ GeV'^ < Sq < (m^ + 0.8)^ GeV'^, the strong coupling constants weakly depend 
on this parameter. 

As an example, let us consider the — )■ S^'^vr^ transition. The dependence of the 
strong coupling constant for the — )■ S^°7r° transition on A'P at different fixed values of 
the So is depicted in Fig. (1). Prom this figure, we see that the strong couphng constant 
for — )■ S^°7r° shows a good stability in the "working region" of M^. This figure also 
depicts that the result of strong coupling constant has weak dependency on the continuum 
threshold in its working region. From this figure, we deduce (?=«()-,() () = 41 ± 7. From the 

b b 

same manner, we analyze all considered strong vertices and obtain the numerical values as 
presented in Table (2). Note that, in this Table, we show only those couphngs which could 
not be obtained by the SU(2) symmetry rotations. The errors in the values of the coupling 
constants presented in the Table include uncertainties coming from the variations of the sq 
and as well as those coming from the other input parameters. Here, we should stress 
that in the Table (2), we only present the modules of the strong coupling constants, since 
the sum rules approach can not predict the sign of the residues of the heavy baryons. Our 
numerical calculations show that the HQET are violated approximately 5% (16%) for the 
coupling constants of the heavy baryons containing 6(c) quark. Finally, we also check the 
SU (3) / symmetry violating effects and see that they change the results maximally about 
8%. 

At the end of this section, it should be mentioned that the predictions of the sum rules 
on light baryon-meson couplings strongly depend on the choice of the structure (for more 
detail see [27]). In connection with this point, here immediately arises the question whether 
or not a similar situation occurs for the case of the heavy baryon-light meson couplings. In 
order to answer this question, we also analyze the coupling constants predicted by the Qn^Jb 
structure, which are presented in Table (2) (see the values in the brackets). From these 
results, it follows that the values of the strong coupling constants of the heavy hadrons 
containing 6(c) quark with light pseudoscalar mesons decrease by a factor of about 8(4) 
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times for the g^ul5 structure. However, sticking on the same criteria as is used in [27], we 
find that the gf_iui)ilb is a more pertinent Dirac structure. 

In conclusion, the strong couphng constants of fight pseudoscalar mesons with heavy 
spin-3/2 baryons have been studied within LCSR. Using symmetry arguments, the Lorenz 
structure independent relations among different correlation functions have been obtained. 
It has been shown that all possible transitions can be written in terms of one universal 
invariant function. Furthermore, it has been observed that the values of the coupling con- 
stants are strongly structure dependent similar to the case of light baryon-meson couplings. 
The numerical values of those strong coupling constants which could not be obtained via 
the SU(2) symmetry rotations have been also presented. 
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Appendix 



In this appendix, we present some details of our calculations, i.e. how we perform the 
Fourier and Borel transformations as well as continuum subtraction. For this aim, let us 
consider the following generic term: 



T= / d^xe'^^ [ due'^'^^f{u) ^''^'^S^L^\ (.1) 
J Jo ^ ' ^ ' 

where Ky is the modified Bessel function of order v appearing in the propagator of heavy 
quark. Using the integral representation of the modified Bessel function 



we get 

T = I d^x C du e^^V(^i)^^^^ii^^ r dt cos(mQt)^= -, (.3) 

where P — p + uq. For further calculations we go to the Euclidean space. Using the 
identity 



1 



T(n) Jo 

we have 



da a"-^e-"^, (.4) 



T = — — — / duf(u) / dt e^'"^* / dy y"^-^ / dv v'^h"'^' / d^xe-'^^-y^''' 

V7rm^r(^) Jo Jo Jo Jo J 



where means vectors in Euclidean space and we will consider only the real part of the 
complex exponential function, e*'"'^*. After performing Gaussian integral over x, we obtain 

^ : " / dufiu) / dte-«*/ dyy^-W .^"V^^f^. (.6) 



0Fm^r( V) Jo ' Jo Jo Jo iv + v) 

The next step is to perform the integration over t. As a result we obtain 



Let us define new variables. 



X^v + y, ^=^^ (-8) 

v + y 
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hence 



T = / duf in) IdXldr A^-V^-(l - r)^-e--I^e-«. (.9) 

Performing Double Borel transformation with respect to the and {p + pY using the 

B{M^)e-'^p' = 5{1/M^ - a), (.10) 

we get 



B{Ml)B{Ml)T = —!^^[\uf{u)[dx[drX'^-'r'^-\l-rr-'e-Ae-- 



m'-gTi^j JO J J 

^ ^ 1 1 - lis 

Performing integrals over u and A, we obtain 

n-\-iy 9 Q 



B{Mt)B{M^)T = / dr/(«o)( — ) r^-^(l - r)^-^e-^^e-^-.% 



''Q^ V 2 



where, Uq = M^+M'i ^ Replacing r = x^, we will have 



(.12) 



"^Q^ l~2~'' •''0 V 4 / 

(.13) 

Finally, after changing the variable r) — and using q'^ — m^, we get 

B(MMM|)T = J.r(^" /(.o)(^J e-^*(^a,/3, Jj, (.14) 



'Q'- y 2 

where 



^2 \ -I POO „2 



M2 y r(Q;) 



with a — and ^ = 1 — i/. 



Now, let us discuss how contribution of the continuum and higher states are subtracted. 
For this aim we consider a generic term of the form 



2 



^ = (M2)«/(«o)^(^«,/3,^J. (.16) 

We should find the spectral density corresponding to this term (see also [28]). The first 
step is to expand /(i^o) as 

f{uo) = T.akul (.17) 
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As a result we get 

[m? + M| j ^''^ i^Mf + Ml ; r(a) A ^"^^ 1) ■ 

Introducing new variables, cri = ^ and cr2 = we have 

k 1 /-oo 

((71 + (72)"+'= r(a) 7i ' / V/ ; 

r{n + k)r{a)Ji Jo 

k poo foo 

r(n + /c)r(a) Jo 

r(n + A;)r(a) Jo V / 



(.19) 

Applying double Borel transformation with respect to Ui ^ — and (72 — >■ — , we obtain the 
spectral density 

X S(s2-{i + rpn'g)). (.20) 

Performing integration over ^, finally we obtain the following expression for the double 
spectral density: 

X e{s,-r)ml), (.21) 



or 

o2 



"(^'•^^^ = ^"- rinV^rw i rf™^-"-'(.-ir'(«:-'<".^)-'-'((s;)'*fe-.)). 

(.22) 

Using this spectral density, the continuum subtracted correlation function in the Borel 
scheme corresponding to the considered term can be written as: 

/ dsi ds2 p{si,S2)e-''/^'e-''/^'. (.23) 

Defining new variables, Si = 2sv and S2 — 2s(l — v), we get 

U'^b = p ds J dv p(si, S2)(4s)e-'^^/^i'e-2^(^-^)/^2'. (.24) 
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Using the expression for the spectral density, one can get 



X 



Safc— r f dv^( (-^)''6(v - 1/2)] 



Integrating over finally we obtain 



- ^'''rin + kna)LJ'2^s^ 



X 



J f2sv/m?Q 



v=l/2 

(.26) 



References 



[1] K. Nakamura et a/., J. Phys. G 37 , 075021 (2010). 

[2] G. Kane, (ed.), A. Pierce, (ed.), "Perspectives on LHC physics", (Michigan U.). 2008. 
337pp. Hackensack, USA: World Scientific (2008) 337 p. 

[3] N. Roberts, M. Pervin, Int. J. Mod. Phys. A 23, 2817 (2008). 

[4] N. Nielsen, F. S. Navarra and S. H. Lee, arXiv: hep-ph/0911.1958. 

[5] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[6] I. I. Balitsky V. M. Braun, and A. V. Kolcsnichcnko, Nucl. Phys. B 312, 509 (1989) 

[7] T. M. Aliev, K. Azizi, M. Savci, Phys. Lett. B 696 (2011) 220. 

[8] T. M. Aliev, K. Azizi, M. Savci, Nucl. Phys. A 852 (2011) 141. 

[9] T. M. Aliev, K. Azizi, M. Savci, Eur. Phys. J. C 71, 1675 (2011). 

[10] T. M. Aliev, K. Azizi, M. Savci, V. S. Zamiralov, Phys. Rev. D 83, 096007 (2011). 

[11] S. L. Zhu and Y. B. Dai, Phys. Lett. B 429, 72 (1998); P. Z. Huang, H. X. Chen, S. 
L. Zhu, Phys. Rev. D 80, 094007 (2009). 

[12] T. M. Ahev, A. Ozpineci, S. B. Yakovlev, V. Zamiralov, Phys. Rev. D 74, 116001 
(2006). 

[13] T. M. Aliev, A. Ozpineci, M. Savci and V. Zamiralov, Phys. Rev. D 80, 016010 (2009). 
[14] T. M. Aliev, K. Azizi, A. Ozpineci and M. Savci, Phys. Rev. D 80, 096003 (2009). 



14 



[15] T. M. Aliev, A. Ozpineci, M. Savci and V. Zamiralov, Phys. Rev. D 81, 056004 (2010). 
[16] T. M. Aliev, K. Azizi, M. Savci, Nucl. Phys. A 847 (2010) 101. 
[17] P. Ball, JHEP 01, 010 (1999). 

[18] P. Ball, V. M. Braun, and A. Lenz, JHEP 065, 004 (2006). 

[19] P. Ball, R. Zwisky, Phys. Rev. D 71, 014015 (2005). 

[20] I. I. Balitsky and V. M. Braun, Nucl. Phys. B 311, 541 (1989). 

[21] K. G. Chetyrkin, A. Khodjamirian, and A. A. Pivovarov, Phys. Lett. B 651, 250 
(2008). 

[22] I. I. Balitsky, V. M. Braun, and A. V. Kolesnichenko, Nucl. Phys. B 312, 509 (1989). 

[23] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Riikl, Phys. Rev. D 51, 6177 
(1995). 

[24] T. M. Ahev, K. Azizi, A. Ozpineci, Phys. Rev. D 79, 056005 (2009). 

[25] A. Khodjamirian, Th. Mannel, N. Offen and Y. M. Wang, Phys. Rev. D 83, 094031 
(2011). 

[26] V. M. Belyaev and B. L. loffe, Sov. Phys. JETP, 57, 716 (1982); Phys. Lett. B 287 
(1992) 176; Phys. Lett. B 306 (1993) 350. 

[27] T. Doi, Y. Kondo and M. Oka, Phys. Rept. 398, 253 (2004). 

[28] V. A. Beihn, A. V. Radyushkin, Nucl. Phys. B 260 (1985) 61. 



15 



